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Abstract 

The hadronic decays r/, rf — ► 3ir and rf — > tjtttt are investigated within the framework 
of U(3) chiral effective field theory in combination with a relativistic coupled-channels 
approach. Final state interactions are included by deriving s- and p-wave interaction 
kernels for meson-meson scattering from the chiral effective Lagrangian and iterating 
them in a Bethe-Salpeter equation. Very good overall agreement with currently available 
data on decay widths and spectral shapes is achieved. 
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1 Introduction 



The hadronic decays of 77 and 77' offer a possibility to study symmetries and symmetry breaking 
patterns in strong interactions. The isospin-violating decays 77,77' — > 3ti, e.g., can only occur 
due to an isospin-breaking quark mass difference m u — m d or electromagnetic effects. While for 
most processes isospin-violation of the strong interactions is masked by electromagnetic effects, 
these corrections are expected to be small for the three pion decays of 77 and 7/ (Sutherland's 
theorem) [1] which has been confirmed in an effective Lagrangian framework [2]. Neglecting 
electromagnetic corrections the decay amplitude is directly proportional to m u — ma- 

Moreover, the 77-77' system offers a testing ground for chiral SU(3) symmetry in QCD and 
the role of both spontaneous and explicit chiral symmetry breaking, the latter one induced by 
the light quark masses. In the absence of 77-77' mixing, 77 would be the pure member 77s of the 
octet of Goldstone bosons which arise due to spontaneous breakdown of chiral symmetry. 

Reactions involving the 77' might also provide insight into gluonic effects through the axial 
U(l) anomaly of QCD. The divergence of the singlet axial-vector current acquires an additional 
term with the gluonic field strength tensor that remains in the chiral limit of vanishing light 
quark masses. This term prevents the pseudoscalar singlet 770 from being a Goldstone boson 
which is phenomenologically manifested in its relatively large mass, m v i = 958 MeV. 

An appropriate theoretical framework to investigate low-energy hadronic physics is provided 
by chiral perturbation theory (ChPT) [3], the effective field theory of QCD. In ChPT Green's 
functions are expanded perturbatively in powers of Goldstone boson masses and small three- 
momenta. However, final state interactions in 77 — > 3ir have been shown to be substantial 
both in a complete one-loop calculation in SU(3) ChPT [4] and using extended Khuri-Treiman 
equations [5]. 

In 77' decays final state interactions are expected to be even more important due to larger 
phase space and the presence of nearby resonances. It is claimed, e.g., that the exchange of 
the scalar resonance a (980) dominates the decays 77' — > tj-kti [6] which has been confirmed 
both in a full one-loop calculation utilizing infrared regularization [7] and in a chiral unitary 
approach [8]. In the latter work, resonances are generated dynamically by iterating the chiral 
effective potentials to infinite order in a BSE, whereas in [7] the effects of the a (980) are hidden 
in a combination of coupling constants of the effective Lagrangian. 

In the present investigation we extend the approach of [8] by including p-wave interactions. 
This will also allow us to obtain more realistic predictions for the decay 77' — > 7r + 7r~7r°, where 
p-waves can — in principle — yield sizable contributions to the decay width and Dalitz slope 
parameters. Furthermore, we study the implications of two very recent experiments by the 
KLOE [9] and the VES [10] Collaborations which have determined the Dalitz plot distributions 
of 77 — > 37r and 7/ — > i]7i + 7i^, respectively, with high statistics. Since these new data have 
not been published yet, we first present the results which we obtain by relying purely on the 
numbers quoted by the Particle Data Group (PDG) [11]. As a second step, we include both 
new experiments separately in the fit and discuss the resulting changes. 

The improved analysis of 77, 77' hadronic decays presented here is also timely in view of the 
planned WASA facility at COSY [12] and MAMI-C [13] which will provide even higher statistics 
for these decays. More precise data will help to constrain the parameters of the Lagrangian and 
pose tighter constraints on the framework employed here. We will illustrate that meson-meson 
scattering phase shifts along with available data on 77, 77' hadronic decays provide a set of tight 
constraints which must be met by theoretical approaches. 
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This work is organized as follows. In the next section details of the effective Lagrangian 
in the U(3) framework are given. Section 3 illustrates our way of incorporating final state 
interactions and includes a discussion of constraints set by unitarity. In Sec. 4 we present 
our results based on data from [11] and the changes which arise if the new, but preliminary 
experimental results by KLOE [9] and VES [10] are included. A critical examination of the 
?7 — ^ 37r data of KLOE based on purely phenomeno logical arguments is presented in Sec. 5. We 
summarize our findings in Sec. 6. 



2 Effective Lagrangian 

In this section we present the effective Lagrangian within the framework of U(3) chiral perturba- 
tion theory and summarize the resulting TT°-r)$-r) mixing [8]. Up to second order in the derivative 
expansion the Lagrangian for the nonet of pseudoscalar mesons n°, K + /~, K°, K°, r] s , rj ) 

reads (note that we do not make use of large- N c counting rules) [14-16] 

£ (o+2) = _ Vo + v^rfQPu) + v 2 (rf X + X ] U) + iV 3 (tf X - X ] U) + V^U* d»U) (tf d^U) , (1) 

where U is a unitary 3x3 matrix which collects the pseudoscalar fields. Its dependence on 7r°, 
rj 8 and rjo is given by 



if 



U = exp | diag(l, -1, 0)"-j + diag(l, 1, "2)^ 



iV2r] 



(2) 



The expression (...) denotes the trace in flavor space, / is the pseudoscalar decay constant in 
the chiral limit and the quark mass matrix M. = diag(m u , ma, m s ) enters in the combination 
X = IBM. with B = — (0| qq |0) / f 2 being the order parameter of spontaneous chiral symmetry 
breaking. The coefficients V* are functions of the singlet field r] and can be expanded in terms 
of this variable 



V, 

v 3 



f 

f 



(0) (2)1$ (4)^ 
u i 1 u i j2 ' i j4 i • • • 



for 



0,1,2,4 



(i)Vo 



-3 y+^3 ? + ... 



(3) 



with expansion coefficients not fixed by chiral symmetry. Parity conservation implies that 
the Vi are all even functions of i] except V 3 , which is odd, and Vi(0) = V 2 (0) = Vi(0) — 314(0) = 
\f 2 gives the correct normalization for the quadratic terms of the mesons. Thus, at a given order 
of the derivative expansion one obtains an infinite string of increasing powers of rjo preceding 
each term in the Lagrangian. 

In order to describe the isospin-violating decays r], rj' — + 3n, we need to work with different 
up- and down-quark masses and their difference m u — m d will be parametrized by means of the 
renormalization scale invariant quantity 



B(m d - m u ) (4) 
which can be expressed in terms of physical meson masses following Dashen's theorem [17] 

(5) 



2 2,2 2 

m K0 - m K± + m w± - m n0 
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up to corrections of 0(e 2 p 2 , {m d — m u )p 2 ). 

While in the isospin limit of equal up- and down-quark masses fr° does not undergo mixing 
with the f]s-Vo system, a non- vanishing quark mass difference m u — m d induces TT -f]s-i]Q mixing. 
The mixing of jr°, 1]$ and 770 is determined by diagonalizing both the kinetic and mass terms 
in the Lagrangian. Since we count the mass of the rj' as a quantity of zeroth chiral order, the 
0(p 4 ) Lagrangian 

£ (4) = I>fa>)C* (6) 

i 

contributes to 773-770 mixing already at leading chiral order [18]. The operators Oi can be found 
in [18, 19], but for completeness we will display in Appendix A those relevant for the present 
calculation. The functions /3j can be expanded in 770 in the same manner as the V* in Eq. (1) 
with expansion coefficients f3if\ At leading order in isospin-breaking and at second chiral order 
the mass eigenstates n°, 77, 77' are related to the original fields fr , rjs, i] by [8] 

7T° = (1 + i?^0 7r o)7T + R^O^T) + RftOrfT)' 
77s = i? 87r 07T° + (1 + R% r ,)ll + R8 V 'V' 

77o = i? 07r o7T° + R 0vV + (1 + R^W (7) 

with the mixing parameters given by 

p (o) _ m e p(0) _ _ p(0) 

„(2) _ p(0) / p(2) 1 2 a \ d(2) _ _p(0) / p(2) . 2 a 

-"-8^0 — -"-8^0 ^-"-TfOTrO -I- 3 ^gmo J , n~ 0r} — n gw0 yrL Sr) -r 3 la 



GMO / ? 



p( 2 ) _ _ -_ _ <]_ pv^; _ _DW i 

0r > ~ ~ V2Pm 2 ' 8r >' " 0r > + V2f 2 

p(2) _ op(0) p(2) p(2) _ 9 p(0)p(2) 

-n- 07r o — ori 87r0 n, 0ri , 7i- 0r? , — Z7x 87r0 7i 8r; ,, 

R (2) tfmHefK + m^) (2) 2^ >17>18 (m; + mJ) 



^K-mJ) (2) _ D (2) S/^m 2 - m 2 ) 



(8) 



7T 7T 



i? (2) 



y 2 " ' Or;' " """" y2 

4/5fm2 + 6/3l 0) (m2+m2) 



/2 

where the superscript on R denotes the chiral order and we have employed the abbreviations 
v? = \P ~ \M\ Pi% = /#> + 3^/2 and ^ >17>18 = 3/?f + /#> - 9/$> + 3/?g>. Note that 
for ^ the mixing parameters given in Eq. (8), which have been derived by diagonalizing 
the Lagrangian, depart from the usually employed orthogonal mixing scheme for 77-77' mixing, 
where i? 8r? / = —Ro v . Isospin-breaking is known to be small, therefore terms of order [m u — m d ) 2 
have been neglected. The quantity = 2v^ / p is the mass of the 77' meson in the chiral 
limit, m 2 = 2_Bt77 and m 2 = 4B(m s + 2m) /3 with 777. = (m u + m d )/2 denote the pseudoscalar 
meson masses at leading order, while the deviation from the Gell-Mann-Okubo mass relation 
for the pseudoscalar mesons is given by 



Am 2 K - ml - 3m 2 6(m 2 - m 2 ) 



i GMO 



m 2 - ml f 2 
where m\ = B{m + m s ). 



P f) _ 6 ^(o) _ 12/3 (°) + 
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Figure 1: Shown is a possible contribution to final state interactions in the decay 77 — > n + n 7r°. 

3 Final state interactions 

One-loop corrections have been shown to be important in the decay 77 — > 37r [4], but even when 
they are taken into account the resulting decay width remains below the measured value [11]. 
In 7]' — > 3tt one expects even larger contributions from final state interactions [8], whereas 
for rj' — * f]7ciT reasonable agreement with experiment can also be achieved in a perturbative 
approach employing infrared regularization [7]. In the present combined analysis of these 
three dominant hadronic decay modes of rj and rf we include final state interactions in a non- 
perturbative fashion as introduced in [8], but extending the work of [8] by taking p- waves into 
account and by improving the fit procedure for the unknown couplings in the chiral Lagrangian 
via Monte Carlo techniques. 

The underlying idea of our approach is that the initial particle, i.e. the rj or //, decays 
into three mesons and that two out of these rescatter (elastically or inelastically) an arbitrary 
number of times, see Fig. 1 for illustration. All occurring vertices are derived from the effective 
Lagrangian and are thus constrained by chiral symmetry. Interactions of the third meson with 
the pair of rescattering mesons are neglected which turns out to be a good approximation, 
particularly for the decays 77 — s- 37r and 77' — > rjim. In the decays under consideration the 
two-particle states either carry one elementary charge or no net charge. Charge conservation 
prevents transitions between these two sets, while the different channels of one set are generally 
coupled. There are nine uncharged combinations of mesons 



For the p-waves there arise some simplifications in the uncharged channels. Due to Bose sym- 
metry contributions from identical particles vanish and the remaining two-particle states can 
be classified according to their behavior under charge conjugation. While for J = 1 7r + 7r~ and 
KK must be C-odd combinations, the other pairs are C-even, so that transitions between the 
two classes of states are forbidden. 

For each partial-wave I unitarity imposes a restriction on the (inverse) T-matrix of meson- 
meson scattering above threshold 



TrV, TT+tT, 777T°, 7)7], K° K° , K + K~ , t/tT , T]'t], rfr)' , 



(10) 



and a set of four charged channels 



7T°7r + , 777T + , K + K°, t/tt 



(11) 



ImT, 



— l 



q, 



(12) 



8nE, 



cm 



5 



with E cm and q cm being the energy and the three-momentum of the particles in the center- 
of-mass (cm.) frame of the channel under consideration, respectively. Hence, the imaginary 
part of Tp 1 is equal to the imaginary piece of the fundamental scalar loop integral G m m above 
threshold ^ 

Gmrnip 2 ) = J ^ _ ~ 2 - ^ _ ^ _ _^ - ^ , (13) 

with m, m denoting the masses of the two particles. In dimensional regularization the finite 
part Gmfn of G mfh is given by 



Gmfh(p 2 ) = Omm(^) + 



16tt 2 



mm m 2 — m 2 m 
1 + In — — H In — 

H A p A m 

■ artanh ■ 



p 2 (m + m) 2 — p 2 

Knrn(p 2 ) = ((m - m) 2 - p 2 ) ((m + m) 2 - p 2 ) , (14) 

where a m m is a subtraction constant which varies with the scale \i introduced in dimensional 
regularization in such a way that G m m is scale-independent [20]. While unitarity constrains 
the imaginary part of the inverse T-matrix, the real part can be linked to ChPT. From the 
effective Lagrangian up to fourth chiral order we derive the c. m. scattering amplitude Afi(6), 
where is the c. m. scattering angle and the indices i, f denote the in- and out-going meson 
pairs, respectively. The decomposition into partial waves reads 

2 

A{6) = ^(2Z + l)A,fi(cos0) . (15) 

1=0 

As usual Pi denotes the / th Legendre polynomial. Note that partial waves with I > 2 do not 
occur, since we consider the effective Lagrangian only up to fourth chiral order. Given a set of n 
coupled channels the partial- wave amplitudes A\ (just as the full amplitude A) are represented 
by symmetric n x n matrices which are functions of the c. m. energy E cm . The symmetry factor 
for two identical particles in a loop is absorbed into the Afi by equipping the matrix elements 
with a factor of 1/V2 for each pair of identical particles in the initial or final state. The inverse 
T-matrix is then given by 

TfW^ + G, (16) 

where the loop integrals G mfh for the different channels are collected in the diagonal matrix G, 
and inversion yields 

T^ll + A-G]- 1 A . (17) 
The expansion of Eq. (17) in powers of Ai ■ G 

T l = A l -A l -G-A l + -- - . (18) 

can then be linked to the loopwise expansion of ChPT. In fact T\ amounts to the summation of 
a bubble chain to all orders in the s-channel, equivalent to solving a Bethe-Salpeter equation 
with A[ as driving term, where all momenta in A[ are set to their on-shell values. 

The partial-wave expansion of the T-matrix may be cast in Lorentz covariant form. For 
the scattering of particles with four- momenta q iy into particles with four-momenta qf, q~f we 
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define the Mandelstam variables s = (qi + qi) 2 = p 2 , t = — g/) 2 and u = (<& — g/) 2 . By means 
of the covariant operators J\ = Ji(s,t — u) 

Js=l, 

t h fj. v i-u (qf - qf)(qj - q}) 
J P = K v q?q f = — + , 

■/d^p 3 , ^yj 

with 

V* = ~9txv + ViiVvIv 2 (20) 
the partial-wave expansions of A and T can be rewritten as 

A = J2iAJi = A s J s + A p J p + A d J d , 
T = J2iTiJi = T S J S + TpJp + T d J d , 

where the Aj, Tj only depend on s. The Ai and TJ are related to the original partial waves Ai 
and T\ by means of a (diagonal) transformation matrix C\ 

Ai = dAiC h Ti = QTiQ, (22) 

and in analogy to Eq. (18) we can write down a Bethe-Salpeter equation for TJ 

Ti = A-TiGU, (23) 

with = C\ GC\. In the present work, we restrict ourselves to s- and p-waves and drop the 
(i-wave part T d . 

We are now in a position to include infinite rescattering processes of meson pairs, which are 
incorporated in T, as final state interactions in the decay amplitudes. In order to introduce 
a common notation for the decays investigated in the present work, we define Mandelstam 
variables 

s = (ph~ Pi) 2 , t = (p h - pj) 2 , u = (p h - p k f (24) 

for the generic process h — > ijk and the p x represent the momenta of the particles. Since 
all decays under consideration happen to have a particle-antiparticle pair in the final state, 
i.e. either 7r + 7r~ or 7r°7r°, which we denote by j and k with j = k, C-invariance dictates that 
the decay amplitude Ahijk(s,t,u) is symmetric under t <-> u. The full amplitude Ahijk, which 
includes s- and p-wave final state interactions, is constructed in such a way that it reproduces 
the tree level result and the unitarity corrections from one-loop ChPT. It reads 

A h ijk(s, t, u) = A hijk (s, t, u) + ( ^ * - u ) [Ti(s) - Ms)] j 

\=o,i ' hi J k 

+ fj2 J i(t^-s)[f l (t)-A l (t)]\ + (j2Ji(u,s-t)[f l (u)-Mu)]] , (25) 

\=0,1 ' hj,ik \=o,l ' hk ^i 

where A hi j k is the complete tree level amplitude from ChPT up to fourth chiral order. The 
differences T\ — A\ are introduced to avoid double-counting of tree graphs and start contributing 
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at the one-loop level. Depending on the subscripts of the parentheses TJ and Ai collect either 
charged or uncharged channels. For identical particles in the final state (7r°7r°) they must 
be multiplied by a combinatorial factor of \/2, in order to cancel the symmetry factor which 
had been absorbed into the potentials. We note that Ahijk does not entail the full one-loop 
result from ChPT, since tadpoles are not included, but they can be absorbed by redefining the 
couplings of the effective Lagrangian [21]. 

3.1 Unitarity 

The approach described above incorporates the relevant pieces of the ChPT one-loop amplitude, 
fulfills unitarity constraints from two-particle scattering and has a perspicuous diagrammatic 
representation of the final state interactions: the summation of a bubble sum in each of the 
three two-particle channels. However, it does not account for three-body interactions in the 
final state, either mediated by the interaction of one of the two rescattering particles with the 
third, spectating particle or by a genuine three-body force. Therefore, the approach does not 
guarantee exact unitarity of the resulting ^-matrix which implies a relation for the imaginary 
part of the decay amplitude Ahijk [22] 

Im A hijk = g ^{^Y^iPi + Pj +Pk-Pa~Pb- Pc) Vbc,ijk Ahabc , (26) 

where T* bci j k represents the complex conjugate of the scattering amplitude for ijk — > abc and 
the sum, which includes the integration over phase space, runs over all possible three-particle 
states which h can decay into. A diagrammatic representation of the unitarity condition is 
shown in Fig. 2. For T a bc,ijk we make an approximation similar to the one already applied to 
Ahijk, i-e., we drop the last diagram on the r.h.s. of Fig. 2 and keep only the graphs involving 
exclusively two-particle rescattering. The first term on the r.h.s. can be expressed in terms of 
the unitarized two-body scattering amplitude T, and Eq. (26) reduces to 

lmA hijk (s,t,u) 

= i E / E { J ^ * - M ') e^*)]* A ^ *> M ') 6(4) ^ +p°- * - 
+i6^E / ^§?E{ J ^- Fr\t)\* A hajc (s, t, u >) + Pc - Pl - Pk )} 

a,c J ^ a ^ c 1=0,1 

E / E { J ^ s> - o [*?**'(«)]• -W s '> u ) 5(4) (pa + Pb - p, - Pj )}, 

(27) 

where, in analogy to the definition of A, we only include two-particle rescattering in the s and 
p partial wave. The Mandelstam variables s', t', u' are defined as 

s' = (Ph - Pa) 2 , ? = (p h - p b f , u' = (p h - p c f . (28) 

The two spectator approximations utilized for A and the r.h.s. of Eq. (27), however, do not 
coincide, since contributions like the one shown in Fig. 3 appear on the r.h.s. of Eq. (27), but 
are missing in A on the l.h.s. The violation of Eq. (27) gives us the possibility to estimate 
the importance of this class of three-body contributions which goes beyond pure two-particle 
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Figure 2: Diagrammatic representation of the unitarity relation in Eq. (26). The crosses indicate 
on-shell particles and phase space integration. 




Figure 3: Contribution which is included on the r.h.s. of Eq. (26), but not in A on the l.h.s. 
The crosses indicate on-shell particles and phase space integration. 
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rescattering as embodied in A. As we will discuss below, it turns out that these deviations 
are rather small for rj —>■ 3tt and rf — > rjim where phase space is narrow and dropping the last 
term in Fig. 2 appears to be a good approximation. Assuming that structures involving more 
complicated iterated two-body interactions or three-body contact terms yield contributions of 
comparable (or smaller) size we conclude that our approach approximates the physical ampli- 
tude reasonably well. For rf — > 3n, on the other hand, where phase space is about a factor of 
seven larger than for rj — > 3tt, neglecting the last term in Fig. 2 is no longer justified and hence 
Eq. (27) is not suited to be an appropriate estimate of unitarity corrections not included in the 
approach. A more detailed study on the importance of three-body effects is beyond the scope 
of the present investigation, but will be addressed in forthcoming work [23]. 

3.2 Isospin decomposition 

In order to study the importance of the various two-particle channels in the final state inter- 
actions and corresponding contributions from well-known resonances like /o(980), a (980) we 
perform a decomposition into isospin channels. Assigning one common mass for all particles of 
an isospin multiplet this can straightforwardly be done for the isospin-conserving decay modes 
rf — > i]7i7i. To this aim, we decompose the interaction kernel of the Bethe-Salpeter equation, 
Eq. (23), into an isospin-conserving and an isospin-breaking part, Aj and A; , respectively, 
so that 

A l = OAf c) T + OAf h) T , (29) 

where O represents the orthogonal transformation which transforms from the isospin to the 
physical basis. Analogous to the definition of TJ we can construct the unitarized amplitude 
in the isospin limit by replacing A\ by A^ in Eq. (23). After substituting in Eq. (25) the 

pieces of the form [TJ — Aj\ by 0[TJ — A\ 1C ^]0 T , the influence of the different isospin channels 
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may be examined by omitting one specific combination of isospin and angular momentum 
quantum numbers. 

The situation is slightly more complicated for the isospin-breaking decays of r\ and rj' into 
three pions. Retaining only one isospin-breaking vertex and inserting it at all possible places 
in the bubble chain, the decay amplitude in terms of isospin channels is found by replacing in 
Eq. (25) the pieces of the form [fj - A{\ by 

0[f/ ic) Gi Af h) Gi f/ ic) - f/ ic) G t Af b) - if b) G t f/ ic) ]0 T • (30) 

While the second and third term in the bracket describe the insertion of the isospin-breaking 
vertex at both ends of the bubble chain, the first one includes insertions at all intermediate 
points. 



4 Results 

We now turn to the discussion of the numerical results of the calculation which are obtained 
from a combined analysis of the decay widths, branching ratios, and slope parameters of the 
considered decays as well as phase shifts in meson-meson scattering. The widths of rj — > 3tt 
and rj' — > i]7i7i have been measured roughly at the 10 % precision level, while for rj' — > 37r° 
the experimental uncertainty is considerably larger and only an upper limit exists for F(t]' — > 
7T + 7r~7r°) [11]. Moreover, some of these decay widths are constrained by the well-measured 
branching ratios 

= T(V - 37r°) = T(rf 3tt°) 

1 T(rj -> tt+tt-tt ) ' 2 T(rj' -> r]7rH°) ' 1 ; 

The Dalitz plot distribution of the decay h — > ij k (with j = k) is conventionally described 
in terms of the two variables 

y/3(u - t) 

x 



2mh{rrih — rrii — 2m 



jk) 

(rrii + 2m jk) \{m h — mA 2 — s] r s 

y = — J 7 LV - 1 , (32) 

2m h m jk {m h - mi - 2m jk ) 

where the m x denote the masses of the respective particles (mj = m k = rrijk) and the Man- 
delstam variables have been defined in Eq. (24). In rj — > 3n measurements (e.g. [24]) a slightly 
simpler definition of y, where (m n o + 2m n +) / m n + is replaced by 3, is usually employed, 

3 \(m h — mi) 2 — s] , s 

y = 9 —\ - 1 ' ( 33 ) 

2m h [m h - mi- 2m jk ) 

but the difference is at the level of 1 % and can be safely neglected. The squared absolute value 
of the amplitude, \Ahijk(x,y)\ 2 , is then expanded for rj' — > r/mi and the charged decay modes 
of r], rf — > 3n as 

\A(x,y)\ 2 =\N\ 2 [l + ay + by 2 + cx 2 + dy 3 + •••], (34) 
while for the decays into three identical particles Bose symmetry dictates the form 

\A(x,y)\ 2 =\N'\ 2 [l+g(y 2 + x 2 ) + ■■■}. (35) 
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For the Dalitz plot parameters a, 6, c of rj — > 7r + 7r~7r° the experimental situation is not without 
controversy. We employ the numbers of [24] , since it is the most recent published measurement 
and the results appear to be consistent with the bulk of the other experiments listed by the 
Particle Data Group [11]. They differ somewhat from the new preliminary results of the KLOE 
Collaboration [9] that has found a non-zero value for the third order parameter d which was not 
included in previous experimental parametrizations. Very recently the Dalitz plot parameters 
of f]' — > i]7i + 7i^ have been determined with high statistics by the VES experiment [10]. While 
we will employ in our fits the experimental Dalitz parameters provided by the Particle Data 
Group [11], we will also discuss further below the modifications of our results when the two new 
and precise (but so far preliminary) data sets [9,10] are taken into account. Note that the slope 
parameters of rj' — > 3n have not yet been determined experimentally, but such a measurement 
is intended at WASAOCOSY [12]. 

From the unitarized partial-wave T-matrix in Eq. (17) one may also derive the phase shifts 
in meson-meson scattering. Hence, our approach is further constrained by the experimental 
phase shifts for -kit — > -kit, KK scattering. The results of the fit are presented in Appendix B. 

The coupled-channels framework entails several parameters, i.e. the low-energy constants 
(LECs) of the chiral effective Lagrangian up to fourth order and the subtraction constants 
a m fn in the loop integrals G mih which are embodied in the coupled-channels T-matrix. It turns 
out that only the fit to the 7T7T phase shifts in the I = J = 1 channel (which includes the 
p resonance) requires a non-zero value of the corresponding subtraction constant a ( k^~ J ~ 1 \fi). 
The regularization scale of G m m is set to fi — 1 GeV for all channels. As a guiding principle 
for the importance of the LECs and in order to reduce their number, we make use of large- N c 
arguments in the effective Lagrangian and set all LECs to zero which are of order 0(1/N%) and 
thus suppressed by at least three powers of 1/N C with respect to the leading coefficients. Their 
effects are expected to be small and can be partially compensated by readjusting the leading 
and subleading coefficients in our fits. Furthermore, we set those parameters to zero by hand 
which turn out to be less sensitive to the processes under consideration. It turns out, that with 
the exception of v ^ and (3$ all parameters of order 0(1/N C ) have a negligible effect when 
varying them within small ranges around zero and can be safely neglected. To summarize, we 
only keep the LECs 

vf\ (3?\ $\ (3?\ $\ p?\ (3[l 0g\ = 0(1) , (36) 
vf\tiS = 0{l/N c ). 

The coefficient is related to the mass of the rj' in the chiral limit, m , and has been con- 
strained to the range 0.00183 GeV 4 . . . 0.00523 GeV 4 in [18], while the rest of the parameters of 
order 0(N l c ) may be varied within small ranges around zero, naturally given by ±N*f 2 / (12A™), 
where n depends on the dimension of the constant under consideration. In conventional ChPT 
the fa term is usually not listed, since it can be absorbed into other contact terms by virtue 
of a Cayley-Hamilton matrix identity. However, this transformation mixes different orders in 
l/N c , hence we prefer to keep this term explicitly, in order to retain the clean large N c behavior 
of the /3j's. 

By fitting to all available (published) data sets of the investigated hadronic rj, rf decays and 
the phase shifts an overall x 2 function is calculated. To this end, we compute x 2 values for all 
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observables i.e. phase shifts, decay widths, branching ratios and Dalitz plot parametrizations, 
divide them by the number of experimental data points and take the sum afterwards. In order to 
find the minima of the overall x 2 function, we perform a random walk in parameter space, where 
only steps which lead to a smaller \ 2 value are allowed, and a very large number of random 
walks with randomized starting points is carried out. We observe four different classes of fits 
which are all in very good agreement with the currently available (published) data on hadronic 
decays, but differ in the description of the decays rj' — > 3n where experimental constraints are 
scarce. 

The errors which we specify in the following for all parameters and observables reflect the 
deviations which arise when we allow for x 2 values which are at most 15 % larger than the 
minimum value. Although this choice is somewhat arbitrary, it illustrates how variation of the 
X 2 function in parameter space affects the results. 



4.1 7/ — ► 3tt 

The results for the decays rj — > 3n, which agree very well with the experimental values, are 
shown in Table 1. Most remarkably our approach is able to reproduce the new, precise value 
of the rj —>■ 37T° Dalitz parameter g measured by the Crystal Ball Collaboration [25] (and 
prevailing the PDG average value [11]) which could not be met in previous investigations [5,8]. 
With regard to [8] this is mainly due to the larger number of chiral parameters taken into 
account in the present work and the improved fitting routine utilized. The PDG number for 
g does, however, not agree with the preliminary g value of the KLOE Collaboration [9] which 
comes close to the results given in [5,8]. The detailed discussion of the new KLOE results for 
the Dalitz plot distributions of rj — > n + n^n and rj — >• 37r° is deferred to Sees. 4.5 and 5. 

When electromagnetic effects are neglected (which is justified according to Sutherland's 
theorem [1]), the isospin- violating decay of 77 into three pions can only take place via a finite 
quark mass difference m u — m d . The decay amplitude is therefore proportional to m 2 e defined 
in Eq. (4) and we have employed the value which follows from Dashen's theorem, Eq. (5) [17]. 
Deviations of the calculated decay widths from the measured numbers could thus be interpreted 
as a hint to non-negligible subleading corrections to the leading order result by Dashen. In order 
to quantify these deviations, one commonly defines the double quark mass ratio 

Q2 = m s -m m s + fa 
m d -m u m d + m u 

and Dashen's theorem yields Qoashcn = 24.1. Differing Q-values lead to decay widths which are 
related to the original one, r Dashen , by 

r= (^^yr Dashcn . (38) 

Taking into account theoretical as well as experimental uncertainties we find from a comparison 
of our results with data Q = 24.0 ±0.6 which is consistent with the result of [8]. Note, however, 
that this obvious agreement with Dashen's theorem merely reflects the fact that our approach 
is capable of reproducing the experimental decay widths of 77 — > 3n. Due to the larger number 
of chiral parameters with increased ranges compared to [8] and the improved fitting procedure 
we can easily compensate the effects from variations in Q by readjusting the chiral parameters 
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T„ o o (eV) 




' 1 




theo. 


422 ± 13 


290 ±8 


1.46 ±0.02 


exp. 


419 ±27 


292 ±21 


1.44 ±0.04 




a 


b 


c 


9 


theo. 


-1.20 ±0.07 


0.28 ±0.05 


0.05 ±0.02 


-0.062 ±0.006 


exp. 


-1.22 ±0.07 


0.22 ±0.11 


0.06 ±0.02 


-0.062 ±0.008 



Table 1: Results for the partial decay widths of 77 — > 37r, the branching ratio 7~i, and the Dalitz 
plot parameters compared to experimental data from [11] and [24]. 

of our approach. We have checked that variations of Q in the range of 20 ... 24 which covers 
the various (and partially contradictory) results in the literature [26] can be accommodated 
within this approach. Therefore, our analysis does not allow for conclusions on the size of the 
violation of Dashen's theorem. 

Extending the work of [8] we have also taken p-wave final state interactions into account. By 
setting these contributions to zero, we find that the decay width of r\ — > 7r + 7r~7r° is reduced by 
a tiny fraction of 0.7% implying rapid convergence of the partial wave expansion. The Dalitz 
plot parameters, which are more sensitive to the precise form of the amplitude than the width, 
are also only moderately altered. Without p-waves we obtain a = — 1.15±0.07, b = 0.29±0.05, 
c = 0.01 ± 0.02. Note that due to Bose symmetry there is no p-wave contribution to the decay 
into three neutral pions. 

Certainly, the most important isospin channel for final state interactions in 77 — > 3ir is the 
1 = s-wave rescattering which is dominated by 7T7T interactions. Omitting this channel reduces 
the decay width by 73%. The other two s-wave channels with isospin one and two, respectively, 
interfere destructively with the former. To be more precise, taking out the 1=1 part, which 
mainly reflects ni] interactions, enlarges the decay width of i] — > 7r + 7r~7r° (77 — > 3n°) by 9% 
(10%), while setting the 1 = 2 channel, which is purely nn rescattering, to zero results in an 
enhancement of the decay widths by 16 % (20 %). The only relevant p-wave contributions arise 
from the 1 = 1 channels irr], irq' which are C-even and thus do not couple to C-odd channels 
related to the p(770) resonance. Neglecting the I = J = 1 channels reduces the 77 — > 7r + 7r~7r° 
decay width by roughly 1 %. The numerical difference to the statement on the importance of p- 
wave contributions in the previous paragraph is due to the fact that for the decomposition into 
isospin channels we use isospin-symmetrized masses, while otherwise we employ the physical 
values of the masses. 

Finally, we have also examined to which extent the amplitude violates three-particle unitar- 
ity in the spectator approximation as described in Sec. 3.1. In order to quantify the violation of 
Eq. (27), we compute the absolute value of the difference between l.h.s. and r.h.s. normalized 
by the modulus of the amplitude Ahijk- Averaged over the whole Dalitz plot we find this vio- 
lation to be (2.5 ± 0.3) % for the process 77 — > 7r + 7r~7r° and — even smaller — (1.3 ± 0.3) % for the 
decay into three neutral pions. The fact that the violation of Eq. (27) is so small is non-trivial 
since only two-body unitarity, but not three-body unitarity is implemented in the definition 
of the decay amplitude. This may suggest that three-body effects (like multiple scattering of 
one particle in the final state with the other two or a genuine three-body interaction) are of 
the same order of magnitude. A more detailed investigation of this issue will be the subject of 
future work [23]. 
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coeff. x 2 y 2 


13.9 ±23.6 


11.8±22.4 


-17.7±9.4 


39.0 ±12.7 


coeff. x 4 


-0.5 ±11.5 


-1.2 ±16.1 


15.4 ±9.8 


-20.5 ±9.4 


poly, order 


6-8 


4-8 


> 8 


> 8 



Table 2: Results for the decay width of rj' — > 7r + 7r~7r° and the leading Dalitz plot parameters. 
The last line denotes the order of the polynomial which is needed to describe the Dalitz plot 
distribution. 



4.2 rj' — ► 3tt 

Only sparse experimental information exists on the decays of rj' into three pions. The experi- 
mental decay width of vj — > 3n° is [11] 

r ( cx P) (?/ _> 3tt°) = (315 ± 78) eV (39) 

which is nicely met within our approach: 

r (theo) ^/ ^ 37r 0) = (330 ± 33) eV _ ( 40 ) 

For the decay into 7r + 7r~7r° only a weak experimental upper limit exists [11] 

r ( cx P) (rj' -> tt+tt-tt ) < 10 keV . (41) 

Due to the large phase space available in these two decay modes of the rj', final state interactions 
are expected to be of greater importance. Indeed we find that in contrast to the processes rj — > 
37r and rj' — > rjivn the Dalitz plot distribution of rj' — > 37r — depending on the choice of the chiral 
parameters — cannot always be well parametrized by a simple second or third order polynomial 
in x and y. Nevertheless, it happens that all our fits may be classified into four groups mainly 
due to the different values of the lower order coefficients in x and y. The numerical results 
for these most relevant coefficients are compiled in Tables 2 and 3 along with the predicted 
width of v( — > 7T + 7r~7r° and the order of the polynomial in x and y which is needed to obtain 
a reasonable approximation to the Dalitz plot distribution resulting from our approach. Note 
that due to charge conjugation invariance, only even powers of x appear. Examples of two very 
different Dalitz plots are shown in Fig. 4. Despite these differing predictions one should keep 
in mind that all fits describe all available experimental data at the same level of accuracy. The 
Dalitz plot distributions of these decays pose therefore tight constraints for our approach and 
must be compared with future experiments at the WASA@COSY facility. 

While in rf — > 37r° p-wave contributions in two-body rescattering are forbidden by Bose 
symmetry, they can be large in rj' — > 7r + 7r~7T° due to large phase space. Interestingly, their 
size varies significantly depending on the cluster of fit parameters. They are largest for the fits 
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cluster 1 


cluster 2 


cluster 3 


cluster 4 


coeff. x 2 , y 2 (V) 


0.1 ±1.7 


-2.7±1.0 


-2.1 ±0.7 


-0.2 ±0.6 


coeff. y 3 


-0.5 ±1.4 


-1.7±0.7 


-0.2 ±0.6 


-0.8 ±0.6 


coeff. x 2 y 


1.6±4.1 


5.0 ±1.9 


0.6 ±1.8 


2.3±1.7 


coeff. y 4 


0.2 ± 1.4 


2.6 ±1.5 


1.6 ±0.8 


-0.1 ±1.1 


coeff. x 2 y 2 


0.4 ±2.9 


5.3 ±2.8 


3.5 ±1.7 


0.2 ±2.5 


coeff. x 4 


0.1 ±1.5 


2.7±1.5 


1.7±0.9 


0.1 ±1.2 


poly, order 


3-6 


5-6 


4-6 


3-5 



Table 3: Results for the leading Dalitz plot parameters of v( — > 37r°. The last line denotes the 
order of the polynomial which is needed to describe the Dalitz plot distribution. 



j I i I i I i L 




(a) (b) 

Figure 4: Sample rf — > tt + 7t~7t° Dalitz plot distribution \A(x, y)/A(0, 0)| 2 of cluster 2 which 
can be described by a fourth order polynomial in x and y (a) and of cluster 3 which must be 
parametrized by a polynomial of eighth order (b). Due to their symmetry under x — > —x only 
the right half of the Dalitz plots is shown. 
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of cluster 4 where setting them to zero diminishes the decay width by 50 % on average. The 
partial width is reduced by 44% (28%) for cluster 3 (cluster 1), while for the parameter sets 
of cluster 2 suppressing the p-wave contributions alters the width by less than 10 %. The large 
higher order coefficients of the rj' — > 7r + 7r~7r° Dalitz plot distribution are mainly due to p-wave 
contributions. If p-waves are omitted, the fits of clusters 3 and 4 can be well parametrized 
by polynomials of fifth order in x and y, while for most fits in clusters 1 and 2 a sixth order 
polynomial would be sufficient, cf. Table 2. Note that in analogy to the decay 77 — > 37T p-wave 
final state interactions with the quantum numbers of the p(770) meson do not occur. 

In i]' — > 3n the contributions from the various isospin channels depend sensitively on the 
cluster, e.g., omitting the I — J — channel in rj' — > 37r° reduces the decay width by 84% 
for cluster 1, while it is enhanced by 132% on average for the fits of cluster 3. For brevity we 
refrain from giving the full list of isospin contributions. 



4.3 T]' — > 1J7T7T 

In Tables 4 and 5 we show the results for the dominant hadronic decay modes of the 77', namely 
the decays into i]7i + 7i^ and rj^n . They are all in very good agreement with the existing 
(and published) experimental data. Furthermore, we have calculated the branching ratio r 2 , 
Eq. (31), which links the two neutral decay modes rf — > 3n° and rf — > rj^n . We find 

r (theo) = (71 ± 7) x 10 _ 4) r (cxp) = ± x ^4 ^ (42) 

and the accordance with experiment is again persuasive. 

In the isospin limit, m u = m d , the decay width r(r/ — > r]7Y + n^) would be exactly given by 
2 T(i]' — > i]7r ir ), due to the symmetry factor for identical particles in the latter process. If, 
however, one is interested in the isosp in-breaking contributions in the amplitude of rj' — > r/mi, 
one ought to disentangle it from phase space effects which are caused by the different masses of 
charged and neutral pions. With an isospin-symmetric decay amplitude, but physical masses 
in the phase space factors, we find a ratio 

r3= r f-^') = 1.78 ±0.02, (43) 

which is smaller than 2 and compares to r 3 = 1.77 ± 0.02 when isospin-breaking is taken into 
account in the amplitude. (For comparison, if the amplitude is set constant and the physical 
pion masses are employed in the phase space integrals, the ratio is given by r 3 = 1.77.) We 
may thus conclude that within our approach isospin-breaking corrections in the rj' — > rjim 
decay amplitude are tiny. The branching ratio r 3 has not been measured directly. If, however, 
we calculate the ratio of fractions Fj/r tota i for these two decay modes using the numbers and 
correlation coefficients published by the Particle Data Group [11], we arrive at 



,( cx p) 



= 2.12 ±0.19 (44) 



by means of standard error propagation. Such a large branching ratio would indicate significant 
isospin-violating contributions in the amplitude. But the experimental uncertainties are sizable 
and should be reduced by the upcoming experiments with WASA at COSY [12] and at MAMI- 
C[13]. 
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^ 7]' —>ri'n-+ w~ 


a 


b 


c 


theo. 


81 ±4 keV 


-0.116 ±0.024 


0.000 ±0.019 


0.016 ±0.035 


exp. 


89 ±11 keV 


-0.16 ±0.06 







Table 4: Results for the partial decay width of rj' — > r\ix+ix and the Dalitz plot parameters 
compared to experimental data from [11]. 







a 


b 


c 


theo. 


46 ± 3 keV 


-0.122 ±0.025 


0.003 ±0.018 


0.019 ±0.039 


exp. 


42 ± 6 keV 


-0.116±0.026 


0.003 ±0.017 


0.00 ±0.03 



Table 5: Results for the partial decay width of rj' — > r]7i TT° and the Dalitz plot parameters 
compared to experimental data from [11]. 

It turns out that p-wave final state interactions are tiny in the processes rj' — > r]7nr. The 
corrections to the decay widths which they generate are smaller than 0.02% and can thus be 
safely neglected. Consequently, in the isospin basis the relevant two-body channels are given 
by s-wave interactions of isospin or 1 states. When examining the influence of these two 
channels on the rj' — > r\ixix partial widths we observe an interesting pattern. By setting the 
1 = channel to zero for the fits of cluster 1 (cluster 2) the widths are lowered by 22 % (22 %), 
while suppressing the 1=1 part reduces them by 81% (72%). When the fit parameters of 
clusters 1 and 2 are employed, the isospin one channel which includes the tail of the ao(980) 
resonance thus appears to be of great importance for the decay mode rj' — > r\ixix confirming 
the findings of [8] and [6]. The situation is, however, reversed if one considers the fits of the 
remaining two clusters. Taking out the 1 = channel in the final state interactions of the 
fits of cluster 3 (cluster 4) diminishes the decay widths by 79% (81%), whereas erasing the 
channel with 1 = 1 reduces it by only 33 % (28 %). Accordingly, for these sets of parameters the 
1 = channel, which incorporates the effects of the /o(980) resonance and the nir correlation 
at lower energies, has higher impact on the decay widths than the a (980) channel. Although 
the fits of all four clusters yield very similar results for all rj' — > rjim observables, the two 
scenarios can be distinguished by their correlation with the processes i]' — > 3-7T provided within 
our approach. Thus, a precise measurement of rj' — » 3ir decay parameters can also help to 
clarify the importance of ao(980) or / (980) resonance contributions to the dominant decay 
mode of the rj' into r/nn. 

The violation of three-particle unitarity as described in Sec. 3.1 is not as tiny as in the 
case of rj — > 37T, but still remarkably small. Using the definition of Sec. 4.1 we find averaged 
deviations of (11 ± 7) % for rj' — > r]ii + 7i^ and (10 ± 6) % for rf — > rjir -!! . It remains to be seen 
whether corrections from other three-body effects which are not included in the approach will 
be of comparable size [23]. 

4.4 Numerical values of the chiral parameters and r\-r\' mixing 

Before presenting numerical results for the chiral parameters we would like to stress that the 
values of the couplings of the effective Lagrangian employed in the coupled-channels approach 
are in general not identical to those in the perturbative framework. First, contributions from 
tadpoles (which include also effects from the so-called on-shell approximation), and t-/w-channel 
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diagrams in the interaction kernel have been absorbed into the coupling constants. Second, the 
BSE summarizes meson-meson scattering in the s-channel to infinite order. The contributions 
beyond a given chiral order are missing in the perturbative approach and lead to changes in the 
values of the couplings when fitting the results to data. Finally, the subtraction point in the 
renormalization procedure can be different in both schemes. Hence, one must expect differences 
in the values of the coupling constants utilized in both frameworks. 

In Table 6 we show the numerical values of the low-energy constants as well as the non-zero 
subtraction constant a^ J ~^ as they come out for the fits of the four different clusters. In 
addition we display the parameters of rj-rj' mixing R 0ri , R Sr) i which are determined by the values 
of the LECs (3f\ and (3[°g in virtue of Eq. (8). Note that compared to the analysis in [8] 
we have increased the number of chiral parameters which — in conjunction with an improved 
fitting procedure — helped to considerably improve the agreement with experimental data on 
hadronic rj, decays. 

According to the mixing parameters the four clusters of fits may be divided into two groups. 
For clusters 1 and 2 Rq v and R 8rj i are both of similar small size and (mainly) positive, while 
clusters 3 and 4 feature a large, positive R 0v and an R 8r] / which is close to zero. Within the 
present analysis the second mixing parameter Rg v >, which characterizes the fraction of the pure 
octet field 77s in the physical rj', turns out to be more tightly constrained by the fit than Rq v 
which describes the singlet content of the r\. In all cases the numerical results for Rq v and 
R 8rj i deviate sizably from an orthogonal mixing scheme, where R 0v = —Rs v '- For comparison, 
a mixing angle of —20° in the one-mixing angle scheme as found in the literature [3] would 
correspond to R 0r) = —R% v ' = 0.34. 

The fitting procedure does not constrain all parameters at the same level of accuracy. While 
some (e.g. 1*3 , (3^\ /3^\ /3^\ /3^) may vary within large ranges (partly compensating each 
other), others like vjf \ (3^\ (3^\ and (3^ are relatively tightly fixed. These boundaries 
constitute important constraints which must be met in future coupled-channels analyses of 
mesonic processes within the approach described here. In particular, the coefficient v^f 1 encodes 
the mass of the 7/ in the chiral limit, m , by virtue of 

2v {2) 

The fact that the 7/ does not become massless in the chiral limit is a consequence of the axial 
U(l) anomaly of QCD which generates in the divergence of the singlet axial-vector current 
an additional, non-vanishing term involving the gluonic field strength tensor. In the effective 
theory this term is represented by . Employing / = 88MeV, the value of the pseudoscalar 
decay constant in the chiral limit [27], we find from the fits of all clusters m = (900 ± 80) MeV 
which is close to the physical mass of the 77'. 



4.5 Recent experimental developments 

Very recently the Dalitz plot distributions of the decays rj ^ 3n and 7/ — > i]7T + 7r~ have been de- 
termined experimentally with high statistics by the KLOE [9] and the VES Collaboration [10], 
respectively. In this section we will discuss the changes of our results when these new and 
precise (though not yet published) data are included in the fit instead of the PDG values. In 
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cluster 1 


cluster 2 


cluster 3 


cluster 4 


t^xlO 3 GeV" 2 


0.82 ±1.65 


0.23 ±1.46 


-1.92 ±0.62 


-1.47±0.83 


^ 2) xl0 3 GeV" 4 


3.15±0.39 


3.21 ±0.49 


3.07 ±0.42 


2.89 ±0.30 


vj 2, xl0 3 GeV" 2 


-0.16 ±0.34 


-0.12 ±0.32 


-0.07±0.17 


-0.13±0.13 




-0.12±0.18 


-0.07 ±0.22 


-0.06±0.19 


-0.02 ±0.32 


^| 0) xl0 3 


-0.47±0.25 


-0.57±0.22 


-0.49 ±0.14 


-0.49 ±0.18 


^ 0) xl0 3 


0.77±0.18 


0.72 ±0.23 


0.69 ±0.19 


0.64 ±0.34 


4 0) xlO 3 


0.05 ±0.55 


0.19±0.51 


0.06 ±0.26 


0.11 ±0.15 


/3| 0) xl0 3 


0.33±0.18 


0.34±0.15 


0.39 ±0.12 


0.42 ±0.14 


4 0) xl0 3 


0.73 ±0.62 


0.86 ±0.66 


0.77 ±0.83 


0.48 ±0.22 


^ 0) xl0 3 


0.00 ±0.30 


0.06 ±0.28 


-0.25 ±0.13 


-0.34±0.15 


$ 0) xl0 3 


0.13 ±0.25 


0.42 ±0.33 


1.01 ±0.49 


0.76 ±0.48 


4 0) xl0 3 


-0.06 ±0.41 


-0.38 ±0.46 


0.02 ±0.51 


0.15 ±0.14 


/^xlO 3 


-0.08 ±0.65 


-0.02 ±0.60 


0.16 ±0.42 


0.23±0.17 


/^xlO 3 


0.08 ±0.35 


-0.03 ±0.31 


-0.21 ±0.22 


-0.25 ±0.14 


/^xlO 3 


0.80 ±0.80 


0.99 ±0.82 


1.51 ±0.47 


1.49 ±0.30 




-6.1 ±0.3 


-6.1 ±0.2 


-6.1 ±0.3 


-6.0 ±0.2 


Ror) 


0.13±0.26 


0.24 ±0.23 


0.61 ±0.13 


0.55±0.17 


R 8v ' 


0.22±0.11 


0.20 ±0.08 


-0.01 ±0.11 


-0.04 ±0.08 



Table 6: Numerical values of the fit parameters itemized according to the four different clusters 
of fits. They also determine the two rj-i]' mixing parameters Rq v and Rs v '- The regularization 
scale in G is set to /i — 1 GeV. 
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rj' — > r]7i + n 




a 


6 


c 


theo. 


—0.116± 0.011 


-0.042 ±0.034 


0.010±0.019 


exp. [10] 


-0.132± 0.019 


-0.108 ±0.033 


-0.046 ±0.022 






a 


b 


c 


theo. 


-0.127 ± 0.009 


-0.049 ±0.036 


0.011 ±0.021 


exp. [11] 


-0.116±0.026 







Table 7: Results for the Dalitz plot parameters of — > ^7T7t when the VES data [10] are 
included in the fit. 

Table 7 we show the results of a fit, where the VES numbers are taken into account. 3 Since 
the amplitudes for rj' — > rjn + n^ and rj' — > 777r°7r° would be equal in the isospin limit and devi- 
ations are thus isospin-breaking and small in our approach, we only include the leading Dalitz 
parameter a of rf — > rjn n° and omit the higher ones which are — assuming only small isospin- 
violating contributions — not quite compatible with the new results of the VES experiment for 
rf — > ?77r + 7r~. Within our approach the c value has a tendency to remain on the positive side 
in contrast to the result of the VES Collaboration, nevertheless our results are in reasonable 
overall agreement with the Dalitz plot parameters extracted from the VES experiment. Most 
remarkably, none of the various other observables (decay widths, branching ratios, Dalitz pa- 
rameters of other decay modes, etc.) is significantly altered when the VES numbers are included 
in the fit, so that the very good agreement of the results with all published data of hadronic rj 
and rf decays is retained. 

Next, we have replaced the Dalitz plot parameters of rj — > n + n^n and rj — > 3n° quoted in 
Sec. 4.1 by the new and precise results of the KLOE Collaboration [9] omitting again the VES 
numbers, in order to avoid interference of these two new, but so far unpublished data sets. The 
results are compiled in Table 8. While it is possible to accommodate the KLOE numbers for the 
a and c coefficients of the rj — > 7r + 7r~7r° Dalitz plot distribution, our results do not agree with 
b and d. In particular, the value of the ^-coefficient b differs from the KLOE number, which 
has been determined very precisely, by more than five standard deviations. Within the given 
boundaries for the low-energy coefficients of the chiral effective Lagrangian our approach is 
unable to produce a b value as small as the number advocated by the KLOE Collaboration [9]. 
Note that such a small value also implies unexpectedly large corrections to the well-known 
current algebra result b = a 2 /4 [4,28]. It may indicate that contributions from higher chiral 
orders of the effective Lagrangian could play a role for this quantity. But the inclusion of such 
higher orders is beyond the scope of the present investigation and will not be discussed here. 

On the other hand, the KLOE result for the leading order coefficient of the rj — > 3tt° 
Dalitz plot, g, cannot be met, while our result still remains compatible with the PDG value, 
— 0.062±0.008 [11]. Generally, we observe the pattern, that reducing the b value correlates with 
an enhancement of the modulus of g. Finally, fitting to the new KLOE numbers destroys the 
agreement of the measured branching ratio r±, Eq. (31), and our result, which is significantly 
increased. The accordance of the rest of the calculated observables with experimental data 

3 Note, however, that the analysis of the VES experiment is still not completed and the quoted numbers may 
slightly change. 
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7] — > 7T + 7T 7T U 




a [9] 


b [9] 


c[9] 


d [9] 


theo. 


-1.054 ±0.025 


0.185±0.015 


0.079 ±0.026 


0.064± 0.012 


exp. 


-1.072 ±0.013 


0.117±0.012 


0.047±0.011 


0.13 ±0.03 




?7 — > 37T 








ri [11] 


theo. 


-0.058 ±0.011 


1.50±0.01 


exp. 


-0.026 ±0.018 


1.44 ±0.04 



Table 8: Results for the Dalitz plot parameters of 77 — > 3ir and the branching ratio r±, Eq. (31), 
when the KLOE data [9] are included in the fit. For simplicity we have added the statistical 
and systematic errors specified in [9] linearly and display symmetrized error bars according 
to the larger value. Note that our coefficients c, d, and g correspond to d, f, and 2a in [9], 
respectively. 

is only marginally affected by including the KLOE results, also the partial decay widths of 
the two i] ^ 3n decay modes which enter r 1 remain consistent with the — admittedly large- 
experimental error bars. 

In Table 8 we employ the r\ value which is determined by the Particle Data Group by 
performing a x 2_ fit using one decay rate and 18 branching ratios (quoted as "our fit" in [11]). 
The result of the most recent direct measurement of r\ [29], however, is a bit larger and has 
also larger error bars: r\ = 1.52 ± 0.12, where we have added statistical and systematic errors 
linearly. Employing this number instead of the PDG value slightly improves the fit to the KLOE 
data, but does not resolve the disagreement with the Dalitz parameters b and g. Taking this 
value for n we find a = -1.049 ± 0.025, b = 0.178 ±0.019, c = 0.079 ±0.028, d = 0.064 ±0.012, 
g = -0.056 ± 0.012, n = 1.51 ± 0.01. 

We mention in passing that after relaxing the naturalness assumption on the size of the 
chiral parameters described at the beginning of this section, we have found a second class of 
fits, which are slightly closer to the results of the KLOE Collaboration for the Dalitz plot of 
77 — > 7T + 7r~7r°. Apart from involving unnaturally large values of some of the LECs they entail a 
g value which is even larger in magnitude than the one of the previous fits, Table 8. Moreover, 
the agreement with the experimental phase shifts of nn scattering in the I = J = channel 
shown in Fig. 5 is considerably worsened. The branching ratio r l5 on the other hand, is not 
altered, cf. Table 8. 

5 Dalitz plot parameters of rj — > 37T 

As pointed out in the previous subsection it is not possible to accommodate the new KLOE 
results for the Dalitz parameters of rj — > 3ir together with the measured branching ratio of 
the two decay modes, r\. In this subsection we will present an explanation how all these 
experimental quantities can be related in a phenomenological way without making use of model- 
dependent assumptions on the construction of the decay amplitudes. 

The main ingredient is the AI = 1 selection rule which relates the rj — > 7r + 7r~7r° decay 
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amplitude A to the amplitude A for rj — > 37r° [4] 

,4(s, i, «) = A(s, t, u) + u, s) + >l(u, s, t) . (46) 

This rule is valid up to tiny corrections from QCD (suppressed by 0(m 2 )) and of electromagnetic 
origin (suppressed by O (a 2 )). 4 In analogy to the experimental parametrization of the Dalitz 
plot distribution, Eq. (34), we assume that the amplitude A can be well approximated by a 
polynomial 

A(x,y) = N[l + ay + f3y 2 + -ix 2 + ■■■] (47) 

with complex coefficients a, f3, 7. We will drop all terms of third order and beyond and work 
with this minimal parametrization of A which is able to describe the experimental Dalitz plot 
distribution as measured by the KLOE Collaboration 5 [9] 

\A(x,y)\ 2 = \N\ 2 [l + ay + by 2 + cx 2 + dy 3 } 
with 6 a = -1.072 ±0.013, b = 0.117 ± 0.012, c = 0.047 ± 0.011, d = 0.13 ±0.03. 

(48) 

Employing the A J = 1 selection rule, Eq. (46), we are able to derive expressions for the 
leading Dalitz plot parameter of 77 — > 37r°, g, and for the branching ratio 77. Since the complex 
normalization factor N is irrelevant for the determination of the g parameter and drops out 
in the branching ratio 77, we are left with six free constants which parametrize the amplitude 
A, the real and imaginary parts of a, (3, and 7. Four of these can be fixed by matching 
\A\ 2 to the central experimental values of a, b, c, d. However, also the remaining two are 
constrained by the fact that the higher order terms x 2 y, x 2 y 2 , x A , and y A , which automatically 
emerge when squaring A, are expected to have small coefficients, since Eq. (48) appears to be 
a good parametrization of the experimental distribution. As an upper limit for the moduli of 
these higher coefficients not observed in experiment we choose the value of the highest order 
experimental coefficient in Eq. (48), d = 0.13. 

Fitting the remaining two parameters in A within the boundaries dictated by the smallness 
of the higher order terms in \A\ 2 to the experimental numbers for 77 and KLOE g we find 

(th«,) = _ .074 ± 0.012, #( cxp ) = -0.026 ± 0.018 [9], 

(49) 

r (thco )= L47±0 .03, r f xp) = 1.44 ±0.04 [11], 

where the theoretical uncertainties represent the propagation of the errors of the input param- 
eters in Eq. (48). While the two numbers for the branching ratio 77 are very well compatible, 
the calculated value of g differs by about two standard deviations from the number extracted 
by the KLOE Collaboration. It is, however, consistent with the value published by the Crystal 
Ball Collaboration [25]. 

We would like to point out that raising the order of the polynomial parametrization of the 
amplitude in Eq. (47) does not alter these conclusions. Although it would increase the number 

4 Our chiral unitary approach iterates isospin-breaking terms and thus includes corrections to the AI = 1 
selection rule, but we have checked that these are numerically tiny. 

5 Note that in contrast to [9] we do not include C-violating terms proportional to x. Therefore, our coefficients 
c and d correspond to d and / in [9], respectively. 

6 For simplicity we have added the statistical and systematic errors specified in [9] linearly and display 
symmetrized error bars according to the larger value. 
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of adjustable parameters, at the same time more and more constraints would be generated by 
the fact that the numerous higher order coefficients of \A\ 2 all have to be close to zero for 
Eq. (48) to be a good parametrization of the experimental distribution. As a matter of fact, 
we have explicitly checked that the inclusion of, e.g., a y 3 term in the parametrization of the 
amplitude, Eq. (47), yields only tiny numerical improvements for the fit to g and r\. As in 
Sec. 4.5, we have verified that our results do not change significantly, when the PDG value 
for T\ is replaced by the most recent direct experimental determination of this branching ratio 
which yields r\ = 1.52 ± 0.12 [29]. Instead of the numbers given in Eq. (49) we then obtain 
g = —0.071 ± 0.012 which is only slightly closer to the KLOE number and T\ = 1.50 ± 0.03. 
The main restriction for the parameters is thus given by the size of the higher order coefficients 
of \A\ 2 and not by the value of r\. 

We have also checked to what extent the phenomenological amplitude described here fulfills 
the unitarity condition discussed in Sec. 3.1. This can be done utilizing purely experimental 
input and thus without making use of unitarized ChPT, since the scattering amplitude TJ 
which enters Eq. (27) may be expressed by the experimentally determined phase shifts of nn 
scattering, see [22] for the explicit expressions. The violation of three-particle unitarity turns 
out to be 10 % (5 %) for 77 — > 7r + 7r~7r° (77 — > 3tt°) on average over the full Dalitz plot when the 
parameters are fixed to physical observables as above, cf. Eqs. (48, 49). The free normalization 
constant N is chosen in such a way that the unitarity violation is minimized at the center of 
the 77 — > 7r + 7T~7r° Dalitz plot. Although the polynomial amplitude does not incorporate any 
constraints from unitarity the violations turn out to be rather modest. If, on the other hand, 
the restrictions on the size of higher order coefficients of |^4| 2 are released and the fit is forced 
to reproduce the central values of the branching ratio r\ and the KLOE g value, we observe 
unitarity violations as large as 43 % for 77 — > 7r + 7r~7r° and 45 % for 77 — > 3n°. 

6 Conclusions 

In the present work we have investigated the hadronic decays 77, 77' — > 3n and rj' — > r/iin within 
a chiral unitary approach based on the Bethe-Salpeter equation. The s- and p-wave interaction 
kernels of the BSE are derived from the U(3) chiral effective Lagrangian up to fourth chiral 
order with the 7/ as an explicit degree of freedom. Within this approach the incoming 77 or 7/ 
decays into three pseudoscalar mesons and then two of these mesons rescatter — elastically or 
inelastically — an arbitrary number of times, while the third meson remains a spectator. The 
final state interaction of the two mesons is described by the solution of the BSE and satisfies 
two-particle unitarity. For the decays 77 — > 3tt and 7/ — > rjim we have also estimated to what 
extend constraints from three-body unitarity, which is not incorporated in the approach, are 
fulfilled and find that the deviations are rather modest. 

The chiral parameters of the approach are fitted by means of an overall x 2 fit to available 
data on the hadronic decay modes of 77 and 77' and meson-meson scattering phase shifts. We 
obtain very good agreement with currently available data on the decay widths and spectral 
shapes. In fact, we observe four different classes of fits which describe these data equally well, 
but differ in their predictions for yet unmeasured quantities such as the 77' — > it + it~it decay 
width (for which there exists only a weak upper limit) and the Dalitz slope parameters of 
77' — > 37T. The results obtained may be tested in future experiments foreseen at WASA@COSY 
and MAMI-C. The hadronic decays considered here along with phase shifts in meson-meson 
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scattering pose therefore tight constraints on the approach and will allow to determine the 
couplings of the effective Lagrangian up to fourth chiral order. It is important to stress that the 
values of the parameters obtained from the fit are in general not the same as in the framework 
of ChPT which can be traced back to the absorption of loops into the coefficients and higher 
order effects not included in the perturbative framework. 

An intriguing feature of the fits is that they accommodate the large negative slope parameter 
g of the decay i] — > 3tt° measured by the Crystal Ball Collaboration [25] which could not be 
met by previous theoretical investigations. This value must, however, be confronted with the 
more recent but yet preliminary g value of the KLOE Collaboration [9]. If we replace the PDG 
data by the KLOE Dalitz parameters of both the charged and neutral 77 — > 37r decay, we do 
not achieve a good overall fit. It appears that the slope parameters of both 77 — ^ 37r decays 
cannot be fitted simultaneously. In addition, fitting to the KLOE data destroys the agreement 
with the experimental branching ratio of both decays which is known to high precision. To this 
end, we have illustrated that utilizing the AI = 1 selection rule which relates both decays and 
taking the KLOE parametrization of the charged decay as input leads in a model-independent 
way to a g value not consistent with the KLOE g result. 

The importance of the various two-particle channels with different isospin and angular mo- 
mentum has been examined as well. For the 77 — > 37r decays we find that the major contribution 
is given by tttt rescattering in the s-wave 1 = channel, while the I = 1,2 channels interfere 
destructively with the former. The p-wave contribution in the charged decay is tiny, since 
available phase space is small and the C-odd channels related to the p(770) resonance do not 
occur. For 7/ — > 7r + 7r~7r°, on the other hand, phase space is considerably larger, and the size of 
the p-wave contributions ranges from 10 % to 50 % depending on the cluster of fits. 

For the decays 7/ — > r\ixix we find that the s-wave 1 = 1 channels dominate for two classes of 
fits which would confirm the importance of the nearby ao(980) resonance as claimed by previous 
investigations. But the other two clusters are dominated by the 1 = channels. These two 
scenarios can be distinguished by their predictions for the 7/ — > 3tt decays. Thus, a precise 
measurement of 7/ — > 37r decay parameters can also help to clarify the importance of ao(980) 
or / (980) resonance contributions to the dominant decay mode of the rf into r/nn. 
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A Fourth-order operators 

For completeness, we tabulate those pieces of the Lagrangian of fourth chiral order which are 
employed in this work. The full list can be found in [18,19]. The fourth-order Lagrangian is of 
the form 

£( 4 )=J>e>;, (A.l) 

i 

where the /3j are functions of the singlet field 770 which can be expanded in terms of this variable 
in the same manner as the V* in Eq. (1) with expansion coefficients (3^\ The operators Oi which 
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are relevant for this work read 



O = (WCCpC,,) , 1 = (C»C,) (CC U ) , 

Q 2 =(C*C v )(C lt C u ), 3 = {CCpCC) , 

O, = - (CCp) (M) , 5 = -(C"C lt M), 

6 = (M) (M) , 7 = (N) (N) , 

8 =\{MM + NN), 13 = -(C»)(C,C»C u ), 



(A.2) 



0u = - (C") (C„) (CC V ) , 18 = - (C") (CpM) , 
where we have made use of the abbreviations 

C^U^d^U, M = U\ + x t U, N = tf X -X ] U. (A.3) 

B Phase shifts of meson-meson scattering 

In Fig. 5 we show the results for the phase shifts of four meson-meson channels: the I = 0, 2 
s-wave and 7 = 1 p-wave mi — > 7r7r scattering as well as 7T7t — > fTA^ with I — J = 0. The 
agreement with the experimental data points is remarkably good. The shaded areas indicate the 
variation of the results when the overall x 2 value of the fits is allowed to exceed the minimum 
by at most 15%, cf. Sec. 4. The variation is particularly small for I = J = 1 scattering 
which entails the p(770) resonance. This fact is, however, not surprising since this channel does 
not enter the hadronic r\ and rj' decays and involves an additional fit parameter, the subtraction 
constant a^' J ~ l \ 
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